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Abstract A convenient and widely used way to study the turbulent plasma in the solar corona is to do statistics of properties 
of events (or structures), associated with flares, that can be found in observations or in numerical simulations. Numerous papers 
have followed such a methodology, using different definitions of an event, but the reasons behind the choice of a particular defi- 
nition (and not another one) is very rarely discussed. We give here a comprehensive set of possible event definitions starting from 
a one-dimensional data set such as a time-series of energy dissipation. Each definition is then applied to a time-series of energy 
dissipation issued from simulations of a shell-model of magnetohydrodynamic turbulence as defined in lGiuliani and Carbond 
(r998), or from a new model of coupled shell-models designed to represent a magnetic loop in the solar corona. We obtain 
distributions of the peak dissipation power, total energy, duration and waiting-time associated to each definition. These distribu- 
tions are then investigated and compared, and the influence of the definition of an event on statistics is discussed. In particular, 
power-law distributions are more likely to appear when using a threshold. The sensitivity of the distributions to the definition of 
an event seems also to be weaker for events found in a highly intermittent time series. Some implications on statistical results 
obtained from observations are discussed. 
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1. Introduction 

The mechanism heating the solar corona to millions of degrees 
remains an open problem, but it is generally understood that a 
great part of the energy dissipation must occur at scales that 
are smaller than the structures that can be resolved by obser- 
vations (x 100 km), as small as 10-100 m (the Kolmogorov 
turbulent cascade dissipation scale). One of the most success- 
ful approaches to reach this four-order-of-magnitude wide gap 
is to assume that the statistics obtained at observable scales are 
still valid at smallest scales. The properties of the global sys- 
tem, from observable to non-observable scales, can then be in- 
vestigated. This is for example the idea underlying Hudson's 
(1991) critical power-law slope of -2 for the distribution of 
flare energies. 

The measurement of the power-law slope for the lowest en- 
ergy flares has in deed been a major goal o f coronal physics 
in the last decade. IXschwanden et al I1I2OOOI) has summarized 
the distributions of event energies that were obtained at wave- 
lengths from X-rays to ultra-violet (UV), and for event ener- 
gies covering a range of eight orders of magnitude from 10'^ J 
("nanofl ares") to 10^^ J ("flares"). It seems — and it is a state- 
ment of lAschwanden et al\ (I2OO0I) — that these distributions 
can be matched together to form a unique power-law distribu- 
tion of slope a; -1.8. 
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However, for the smallest events, mainly observed in UV 
by filter imaging instruments, some observations seem not 
to mutually agree, and they seem not to fit in the global 
distribution. A possible expla nation for this is given by 
Aschwanden and Charbonneaul (I2OO2I) : the energy of an event 
could be wrongly deduced from the observable quantities (like 
intensity in some spectral lines), leading to a systematic error 
in the distribution of event energies. Another explanation could 
be that all authors do not agree on what they mean by "event", 
i.e. the fact that inequivalent definitions exist in the literature. 

Indeed, defining an event is likely to be much more dif- 
ficult for low energy events than for high energy events. For 
high-energy events, whose distributions are in general derived 
from X-rays and radio observations, there is little ambiguity on 
what is an event: events are very rare (10^^ s"' for the whole 
solar disk between 10^^ and 10^"* J) and well-separated by long 
low-flux times. On the contrary, low-energy events can be very 
close in space and time, making it difficult to separate them, 
either because they occur on the same line-of-sight or because 
they are smaller than the instrumental resolution or shorter than 
the time resolution. The difficulty is even bigger if we sub- 
scribe to the idea that the corona is in a self-organized criti- 
cal state so that small events trigger other events, leading to 
avalanches as illustrated by th e sandpile paradigm dBak et al\ 
I1988HLu and H amilton" 1991*): among all these events, which 
ones should be used to do statistics? 
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We think that this difficulty has been underestimated when 
statistics obtained from observations or simulations by different 
methods have been compared. In fact, the definition of an event 
which is used is very rarely discussed (contrary to the influence 
of the relationship between the observable parameters and the 
physical variables of an event), and is sometimes even not given 
precisely. 

In this paper we give some definitions that could be used, 
mainly inspired from definitions which have already been used 
in the past. We choose to restrain ourselves to events defined 
from a one-dimensional data set, namely a time series of en- 
ergy dissipation, so that the definitions can be easily com- 
pared. We then produce statistics of events (histograms of 
event energies, durations and waiting times), for different def- 
initions of an event, and compare them. The time series we 
use are the data output by a shell-model of MHD turbulence 
jGiu liani and Carbone, I l998l) . and by a new model of coupled 
shell-models describing Alfven turbulence in a coronal loop 
( iBuchlin et aL.200 4). However, the aim of this paper is not to 
study shell-models of MHD, but to see to which extent the defi- 
nition of events influences their statistics, even in a simple case 
of events detected in a time series of energy dissipation. 

2. Event definitions 

We present here a basic list of possible definitions of an event 
when a "signal" e(f), which is the time series of the power dis- 
sipated in the system, is given (Fig.[0. Most of the ideas of this 
list come from the definitions that have been chosen in papers 
found in the literature. For each event, we get: 

- E, the total energy dissipated during the event, 

- P, the peak power of energy dissipation, 

- T, the duration of the event, 

- te, the time of the event, necessary to obtain the waiting 
times T„„ i.e. the (quiescent) time between two consecutive 
events. 

2.1. Peaks 

Definition 1 (peak). An event corresponds to a local maxi- 
mum e(f,„) in the signal e(f). The time of the event is = t„„ 
the peak dissipation power is P = £(te), the total dissipated en- 
ergy is £ = ' e(f) df where and are the two local minima 
around t^, and the event duration is T - ti, - ta. 

Variant 1.1 (peak-background). The background b(t) is the 
affine function defined between the points (fa,e(fa)) and 
(tb, e(tb)). With the notations of definition ^ the time of the 
event is t^, the peak dissipation power is f - b{tg), the total dis- 
sipated energy is £ - b(t) dt - E - (eitg) + e(ti,)) ■ T/2, and 
the event duration is T. 

2.2. Threshold 

Definition 2 (threshold). A threshold ethi is chosen, and an 
event is a part of the signal e(f) which stays above ethr: more 
precisely, it is a maximal connex part V = [ta, tb] of the set 



[t I e(f) > ethi). The total dissipated energy is E - e(t)dt, the 
peak dissipation power is f = maxy e(t) and the event duration 
is r = tb - ta- The time of the event is the time at which the 
maximum of €{t) on V is attained: e(fp) = maxy e{t). 

There are several other alternatives to define the time of 
the event, like te - {tb + ta)/2 (the middle of interval V) or 
te - ^ Jy e(t) ■ t dt (the barycenter of the event, weighted by e). 
But these variants do not change the statistics of P, E, and T, 
and they seem to have little influence on the statistics of t^,. 

We can consider the threshold level ethr as a background 
level, giving the following definition: 

Variant 2.1 (threshold-background). Using definition |2l and 
its notations, the time of the event is f^, the peak dissipation 
power is P - ethr, the total dissipated energy is £ - ethr dt = 
E - ethr ■ T, and the event duration is T. 

The threshold can be chosen as a function of the overall av- 
erage e and standard deviation o-g of e{t). It can also be chosen 
iteratively, by using the average and standard deviation of the 
time series during the quiescent times between events (which 
h ave been defined by th e previous iteration of this process), as 
in lBoffettaef "011(119991) . 

2.3. Wavelet analysis 

This method assumes that we have built the time-scale plane 
^(fo, s) for e(f), by convolution of e(f) by the wavelets w,„,j(?) = 
l/s ■ wo{{t - f())/i). A mother wavelet wq which seems well 
adapted to the shape of events is the second deri vative of a 
Gaussian ("Mexican hat"). When the noise is in 1 / f . ISanz et al.\ 
( ^001) have shown that the Mexican hat is the best wavelet to 
find enhancements of the signal. 

Definition 3 (wavelet). An event corresponds to a local maxi- 
mum y(te, Se) in the time-scale plane ^(fo, s). The time of the 
event is fg, its duration T is the scale Sg, its total energy is 
E = y(te, Se). Its peak power P can be defined as maxy e(f) 

with V - [te - Sell, te + Sejl]. 

For better accuracy in the definition of Se and te, we need to 
have a good resolution in the time-frequency plane, i.e. we use 
a continuous wavelet transform. As a result, and also because 
we have used the Mexican hat wavelet, the wavelets are non- 
orthogonal, and there is a risk to find several non-independent 
events where there is only one event. To avoid this, it is pos- 
sible to impose that two events are separated by at least some 
distance in the time-scale plane (typically a factor 2 in scale 
and an interval Se in time), but in practice this is not necessary 
as the time-scale plane is sufficiently smooth. 

3. IVIodel time series 

3.1. Shell-model 

The results of the next section are based on data output from a 
shell-model of incompressible MHD turbulence. In such mod- 
els (iGloaguen et aL.1985:.Biskama.l994:.Frick and SokoloS 
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Figure 1: Definitions of events and of event characteristics (t^, T, P, and E). Left: peaks (definition [Q, with (top) and without 
(bottom) background detection. Middle: threshold (definition with (top) or without (bottom) taking the background into 
account. Right: peaks in the wavelets time-frequency plane (definition |3}- For this definition, events are marked as crosses in 
the time-frequency plane, indicating the event time and duration T. The total energy E is the wavelet amplitude (color) at this 
position. 



Il998t iGiuliani and Carbonel Il998l) . the Fourier wavenumber 
space is divided in concentric shells 5„ = {A: | \\k\\ e [kn,k„+i]} 
with k„ - koA" and usually A - 2. A single complex scalar 
value u„ represents the velocity increments \u(x + {) - u{x)\ 
on scales { for Injt e S„. The same holds for the scalar 
value bn representing the magnetic field increments on the 
same scales (. This model is the magnetohyd rodynam i c ana- 
log of the GOY (Gledzer-Ohkitani-Yamada: iGledzeil Il973t 
LYamada and Ohkitani, 1987, 1988a^b) shell-model of fluid tur- 
bulence. It is governed by the following equations, given in 
iGiuUani and Carbonel ( 1 19981) : 

dZ' 

= -kl{v^r„ + y-Z„^) + ik,j: + /„' (1) 

where Z„' - v„ + sb„ are the Elsasser variables, 5 = ±, 
V* = (v H- sT])/2 are combinations of kinematic viscosity and 
resistivity, are external driving forces, and is the term cor- 
responding to local non-linear interactions between shells. For 
a given shell n, this term involves the neighbors and second- 
nearest neighbors of the shell n, modelling local interactions 
between triads of consec utive modes. The detailed coe fficients 
of this term are given in iGiuliani and Carbonel lll998l) and de- 
pend on the dimensionality {e.g. 2D or 3D) of the physical 
MHD system the shell-model represents, via the conservation 
of the MHD invariants. 

This model can describe the evolution of modes over a 
wide range of wavenumbers with just a few dozens of de- 
grees of freedom. It is thus very interesting for studying 
MHD turbulence with high Reynolds numbers, and intermit- 
tency. It actually exhibits typical properties of MHD turbu- 
lence, from wide power-law spectra to dynamo effect in 3D 
tGiuliani and Carbone. .1998.) . including spatial and tempo- 



ral iri termittency ( IGiuliani and Carbon^ Il998t iBoffetta et all 
l2002h . 

The equations^are solved numerically and we get the time 
series of dissipated power e(t) - Xi„ A;^(v|M„(f)P + /7|fe„(0P), 
which is our variable of interest. To obtain the first time se- 
ries shown on Fig.|2 hereafter known as the time series <1>, 
we used 24 shells (representing k - I to k x 8.4 10^), with 
A - 2, V = 77 = 10"". We performed 10^ variable timesteps 
(determined by a CFL condition) with a 3"^-order Runge-Kutta 
numerical scheme. The time series <2>, also shown on Fig.|2| 
was obtained with the same parameters, except that the dissipa- 
tion coefficients v and t] where ten times higher than for <1>. 

3.2. Coupled shell-models 

Section |5] uses also data from a version of a shell-model de- 
signed to model a region of space where a dominant magnetic 
field B() exists, like in a coronal loop (Buchlin et al., 2004). In 
this model, shell-models of 2D MHD are coupled by Alfven 
waves travelling along Bo, and energy is only input by move- 
ments of the photospheric footpoints of the loop. This geomet- 
ric setup is t he same as the one use d for the cellular automaton 
described by 'Buchlin et al. ' ('2003'), and it gives a model simi- 
lar to the one described by Nigro et al. (2004). Here we use an 
independent implementation of these ideas to obtain the time 
series <3>. 

3.3. Characteristics and iritermittency of the time 
series 

All these time series were rescaled so that their average e is 1 
and are shown on Fig|2 Their basic characteristics are summa- 
rized in table[2 From time series <1> to <3>, the ratio of the 
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Figure2: From top to bottom, time series <1>, <2>, and <3>. 



Table 1: Summary of the characteristics of all 3 time series, 
which were normalized so that their average is e = 1 : number 
of data points, number of peaks, standard deviation, and maxi- 
mum value. 





Data points 


Peaks 




^max 


<1> 


453,628 


51,507 


1.98 


169 


<2> 


985,162 


56,136 


2.45 


420 


<3> 


1,000,000 


305,738 


11.33 


1971 



maximum (or the standard deviation) to the average grows, and 
longer, quiet times exist between the intervals with higher dis- 
sipation. It seems that intermittency is higher for <3> than for 
<2>, and that it is also higher for <2> than for <1>. This is ver- 
ified by plotting the flatness' of these time series as a function 
of the temporal scale (Fig.|3}. 



4. Comparison between statistics for different 
event definitions 

In this section, starting from the time series e(f) number <1> 
produced as explained in section ITTI we compare the effect of 
the definition of an event for the following statistics: 

- normalized histograms (i.e. experimental Probability 
Distribution Functions - or PDFs) of event peak dissipation 
power, total energy and duration, as defined by the difi'erent 
event definitions we use, 

- PDFs of waiting times between events, i.e. the time be- 
tween two successive events. This corresponds to the lami- 
nar, quiet time between events. 



' We use the following definition for the flatness F(t) of the time 
series 6(f): F(t) = 5^(t)/(5^(t))-, where 5'(t) = Mt + r) - 6(01"), is 
the structure function of index q for the timeseries e. 
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Figure 3: The flatness of the three time series increases when 
the scale 6t decreases. This behavior is a signature of intermit- 
tency, and is stronger for the time series <3> and lower for time 
series < 1 >. Note that the flatness remains much higher than the 
Gaussian value 3 even at large time scales, as a consequence of 
the non-Gaussian distribution function of the values taken by 
the time series. 



4.1. Peaks 

As each peak of the time series is counted as an event, defi- 
nitions ^ and ^3 give a lot of events, even in the case of our 
numerical data, which has no noise: for time series <1> for ex- 
ample, one data point over nine is a local maximum, and corre- 
sponds thus to an event. When noise is present, a smoothening 
of the data at the scale of the shortest events may be necessary 
before searching for events. Furthermore, the set of the events 
is a partition of the time series (the end ti, of one event is the 
beginning of the next event), all the energy of the time series 
is contained in the events: £, - J e(t) dt. 

The distributions of P, E, T and t^. (Fig.0} have approxi- 
mately the same shape, which is neither a power-law, nor an ex- 
ponential or gaussian. The tail of the waiting-time distribution 
could even be fitted to an exponential (Fig.|5}, in contradiction 
with previous studi es of shell-models, which used another de f- 
inition of an event teoff^etta et aZillQQgtlLepreti et a/.Ll200ll) . 

With these definitions, even the smallest peaks are counted 
as events, and this breaks the waiting times into small parts, 
leading to a cut-off' of the tail of the waiting-time distribution. 
To decrease this effect, we may exclude the smallest events (e.g. 
those with a peak power lower than a given threshold) from the 
analysis, which gives the following variant of definition[3 

Variant 1.2 (peak-threshold). A threshold eth, is chosen. 
An event corresponds to a local maximum e(t,n) in the signal 
e(f), provided that e(t„,) > eihr- The time of the event is te - tm, 
the peak dissipation power is P = e{te), the total dissipated en- 
ergy is £ = e(f) dt where ta and th are the two local minima 
around f^, and the event duration isT = tb - ta- 

Note that this is not the same as using definition |2] for a 
given threshold ethi, on a maximum connex part V of {f | e(f) > 
ethr), definition 1 .2 will find as many events as there are peaks 
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Figure 4: Events statistics of time series <1> for definition 
(peaks): peak power of energy dissipation, dissipated energy, 
duration, and waiting times. The straight dashed Hne corre- 
sponds to one event per histogram bar; as the histogram bars 
are spaced exponentially, its slope is - 1 . The dashed histograms 
are an estimation of the discretization error when building the 
histogram, computed assuming Poisson statistics in each bar 
(+ "^V where is the number of events in a given histogram 
bar). 




C.I r 



0.00 0.02 0.0 8 
Wolting time 

Figure 5: Waiting-time distribution of events from time series 
<1> for definition[2(peaks), in lin-log scale. 



of e(f) on the interval V, whereas definition|2]will find only one 
event. 

As a result (Fig.|6j, it is clear that the PDF of P is cut below 
the value of ethr, with no modification of its shape: this means 
that only the tail f > ethr = 1 of the histogram of P in Fig. |3 
is left, and this tail could be fitted to a quite narrow power-law 
of slope -2.95. The PDFs of E and D do not change dramati- 
cally, except that the left part is weaker because of the correla- 
tions between P, E and D. The most interesting effect of using 
variant 1 .2 instead of definition^is for the waiting-time distri- 
bution: it exhibits now a clear power-law of index - 1 .98 over 
2.5 decades. This is made possible by the fact that small wait- 
ing times associated to small events in the case of definition[2 
are now replaced by a smaller number of long waiting times, 
leading to a reinforcement of the right part of the histogram of 

Tw 



Figure 6: Events statistics of time series <1> for definition 1.2 
(peak-threshold) with ethr = 1 (average of time series). The 
statistics for definition ^ (peaks) are shown as a reference, in 
dashed lines. 

4.2. Threshold 

With definition |2] histograms of P, E, T, and are quite 
clearly power-laws (Fig.^}, even if they are not very wide for 
this weakly intermittent time series. The slopes of these power- 
laws are 2.89 + 0.06 for P, 2.3 1 + 0.05 for E, 2.48 + 0.06 for T, 
and 1.79 + 0.02 for T„. 

These power-law tails still exist when the threshold is con- 
sidered as a background and is removed (variant ITTI Fig.|8}, 
but the left part of the histograms is then almost flat in log- 
arithmic axes. This is for example quite straight-forward for 
the distribution of P, as removing the background is shifting 
— in linear axes — the distribution of P to the left. However, 
the right tail of the distributions, e.g. for P » ethr remains al- 
most the same when the background is removed. It seems that 
background removal does not help understand the statistics of 
events. 

Methods using a threshold are very widely used when 
events are searche d in time series, as w ell from numeri - 
cal simulations: .Dmitruk et al. (1998); .Einaudi et al. 1 (Il996h : 
iGeorgoulis^f «/■! lll998ll (2D RMHDV [RoffenB et al.\ (Il999l) 
(MHD shell-models) as from X-rays observations: Pearce e t al\ 
(1993); Cix )_sby_ef aL (1993); Wheatland ela^- (1998). It actu- 
ally seems to be well-adapted to instrumental constraints of 
sensitivity and noise levels. 

The drawbacks of this definition are that it misses the 
lowest-energy events (leading to a cut-off of the left part of 
the energy histograms), and that it cannot separate close high- 
energy events. This definition is also not adapted to non- 
stationary time series: in this case, the threshold should adapt 
to the local statistical characteristics of the time series. 

4.3. Wavelet analysis 

Definition|3]produces the histograms of P, E, D, and t,v shown 
on Fig. |9] The histogram of events durations is a power-law 
over more than 2.5 decades. The wide and flat left part of the 
histograms of P and E, which include events much smaller than 
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Figure?: Events statistics of time series <1> for definition |2] Figure 9: Events statistics of time series <1> for definition |3] 
(threshold). (maxima in wavelet time-scale space). 
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Figure 8: Events statistics of time series <1> for definition 12. II 
(threshold-background). Only the distributions which are dif- 
ferent from the one obtained by definition|2l(threshold, Fig.Q 
are represented. 

with other definitions, suggest to use a variant of definition |3] 
similar to variant 1.2 of definition[2 where the smallest events 
are simply not taken into account: 

Variant 3.1 (wavelet-threshold). A threshold isthr is chosen. 
An event coiTesponds to a local maximum y(te, s^) in the time- 
scale plane y(to, s), provided that y{te, Se) > E,hr. The time of 
the event is te, its duration T is the scale Se, its total energy is 
E = y{te, Se). Its peak power P can be defined as maxy e(f) with 

V^[te-Se/2, te + sJ2]. 

As for definition 1.2 (peak-threshold), the distributions of 
P, E and D do not change much, but a power-law is recovered 
(Fig.llO> for the waiting-time distribution. 

5. Intermittency and sensitivity to event definition 

5.1. Sensitivity to event definition 

Now we use all the 3 time series described in section|3l It seems 
that the distributions of event energies E obtained by definitions 
n]and|3lare closer from the power-law obtained by definition|2l 
in the case of time series <3> (Fig. II It) than in the case of time 
series <1> (Fig. II Ih). The waiting-times distributions (Fig.ll2> 
display the same behavior. In general, distributions obtained 
from higher-intermittency time series seem to be less sensitive 
to the definition of an event than low-intermittency time series. 
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Figure 10: Events statistics of time series <1> fordefinition l3.1l 
(wavelet-threshold) and isthr - 1 ■ The statistics for definition|3] 
(maxima in wavelet time-scale space) are shown as a reference, 
in dashed hnes. 

5.2. Sensitivity to tliresliold, for definition\^(tiiresiiold) 

In the case of events defined by a threshold (like definition 
13, the slope of event energy histograms may depend on the 
threshold. Here we choose different values of ethr between 
and € + 5o-^ where e is the time series average (e = 1) and cr^ 
is the standard deviation shown in table[2for each of the time 
series. As a result, the number of events (Fig. 113b ) is 1 when 
the threshold is ethr = (the whole time series is one event); 
it increases to a maximum, attained between e and e + cr^ , de- 
pending on the time series characteristics; then it decreases (ul- 
timately, the number of events is when ethr > fmax, where emax 
is the maximum value of the time series). 

Figure ^] shows the power-law slope of the histograms 
of P, E, T, and T^^. as a function of the normalized threshold 
(fthi - e)/o"e- In general, time series <1> and <2>, which come 
from the same simple shell-model and which are less intermit- 
tent than the time series <3>, follow quite the same path, (a) 
The distributions of peak dissipation power P have a slope ^ 2 
for a low threshold ethr, and become steeper when ethr increases. 
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Figure 11: (a) Total energy distributions for events determined 
by definitions n (peaks; dotted line), |2] (threshold; plain line), 
and |3] (maxima in wavelet time-scale space; dashed line), for 
time series <1>, (b) Same figure for time series <2>, (c) Same 
figure for time series <3>, All plots have the same scale. 
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(b) <2> 





Figure 13: Slope of the peak power (a), total energy (b), du- 
ration (c), and waiting time (d) distributions, for events de- 
fined by definition |2] (threshold), as a function of the normal- 
ized threshold (ethr - e)/o"e- The plain, dashed, and dotted lines 
correspond to time series <1>, <2>, and <3> respectively, (e) 
Number of events detected as a function of the normalized 
threshold, (f) Proportion of the duration contained in events, 
as a function of the normalized threshold, 

fact that the deviations of time series <3> are larger than for 
the other time series. 



'_ ■ ^ 
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Figure 12: (a) Waiting-time distributions for events determined 
by definitions [2 (peaks; dotted line), |2] (threshold; plain line), 
and |3] (maxima in wavelet time-scale space; dashed line), for 
time series <1>, (b) Same figure for time series <2>, (c) Same 
figure for time series <3>, All plots have the same scale. 

The slope for time series <3> is slightly more sensitive to ethr 
than the other time series, (b) The slope of the distributions of 
energy E also increase with ethi, except for time series <3>, 
for which it is almost constant, (c) The statistics of the dura- 
tions T exhibit the same features than the statistics of E. (d) 
On the contrary, the slope of the distributions of the waiting 
times T„, decreases when ethr increases. Again, it is almost con- 
stant for time series <3>, (f) The proportion of the time se- 
ries duration contained in events decreases when the threshold 
increases. This decrease is stronger for the lowly intermittent 
time series. 

Time series <3> seems to be the least sensitive to the value 
of ethr- Note that, by using thresholds expressed as a function 
of instead of absolute thresholds, we have taken care of the 



6. Discussion 

We have investigated the dependence on the definition of 
"events" of the statistics of events obtained from an energy dis- 
sipation time series. Not very surprisingly, the statistics of peak 
power, energy content, duration and waiting times of events dif- 
fer when different definitions are used. 

Especially for low-intermittency time series and for 
waiting-time distributions, power-law distributions are recov- 
ered only when a threshold is used, either when searching for 
events (definition^ or after having searched for events by an- 
other means (definitions 1,2 and 13, U , It is also interesting to 
note that the waiting-time distribution, whi ch is used to test 
the Poissonian nature of the flarin g process (Wheatland et M 
J398: Lepreti et ai. 2001: Wheatland and Litvinenko, 2003) . 
can have a power-law or an exponential tail, depending on the 
definition of events. 

For observational studies, where the smallest events are av- 
eraged over the line of sight and the spatial and temporal steps, 
some of the intermittency is lost. In this case we need to use 
a definition which gives statistics as close as possible from the 
statistics of the underlying (non-averaged) signal (which is in- 
termittent enough for events statistics to be almost indepen- 
dent from their definition). The definitions using a threshold 
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seem to be adequate from this point of view. The presence of 
noise in observations gives also a strong support to this kind 
of definitions. However, these definitions have also drawbacks, 
in particular the difficulty of choosing a threshold for a non- 
stationary time series. 

Other definitions like |3 which uses wavelets can have in- 
teresting properties separating simultaneous events at different 
scales, but the smallest events obtained by this means seem to 
be not significant. Alternatively, events could be defined iter- 
atively from the time-scale plane: the first event is defined by 
the overall maximum of the time-scale plane, the correspond- 
ing wavelet is subtracted from e(f), a new time-scale plane is 
computed, and this process is done again to find each of the 
next events. Local Intermittency Measure (LIM: Farge, 1990) 
could perhaps also be used for this purpose. However, these 
ideas have not been investigated further yet and an iterative def- 
inition may be computationaly very expensive compared to the 
other definitions. 

Let us now return to the motivation behind the determina- 
tion and discussion of event properties and statistics for coro- 
nal physics. There are two main reasons for these studies, es- 
sentially related to bridging the gap between observable time 
and spatial scales and the sub-resolution physics. On the one 
hand, we would like to understand whether analogous physical 
processes, namely flares, conserve scale-invariant properties at 
unobservable scales and are responsible for the existence of the 
quiet corona as we know it. On the other, one would like to link, 
as far as possible, large scale physical models and numerical 
simulations to the observations without reproducing in detail 
the microcosm of a single small-scale event (though this may 
be desirable and necessary for the largest scale manifestations, 
such as for example the Bastille day flare), but by comparing 
global statistical properties. 

When searching for the answer to the first question, one 
must clearly use an event definition which conserves the to- 
tal energy in the signal, as one is searching for a quantitative 
confirmation (again, much care is needed, since the average 
corona exists to some extent precisely because we are unca- 
pable of observing fluctuations at sufliciently small energy and 
time-scale, i.e. it is by definition a background). In the second 
case however, where there are undoubtedly large differences 
between numerical models and observations in the richness of 
the physics and dynamical range, one must be careful to ana- 
lyze events in the way most appropriate to glean characteristic 
properties of the fluctuations and turbulence at the available 
scales. Hence the requirements are to go beyond simply the 
energy distributions of events and analyze other characteristic 
features such as anisotropy of the spectra, intermittency, and 
higher order structure functions of the fields. 
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